arXiv: 1509.07110v2 [cond-mat.supr-con] 18 Mar 2016 


Bounds on Nanoscale Nematicity in Single-Layer FeSe/SrTiOa 


Dennis Huang,^ Tatiana A. Webb,^’^ Shiang Fang,^ Can-Li Song,^’* 

Cui-Zu Chang,^ Jagadeesh S. Moodera,^’^ Efthimios Kaxiras,^’^ and Jennifer E. Hoffman^’ 

^Department of Physies, Harvard University, Cambridge, Massaehusetts 02138, USA 
^Department of Physies & Astronomy, University of British Columbia, Vaneouver, British Columbia V6T IZl, Canada 
^Franeis Bitter Magnet Laboratory, Massaehusetts Institute of Teehnology, Cambridge, Massaehusetts 02139, USA 
^Department of Physies, Massaehusetts Institute of Teehnology, Cambridge, Massaehusetts 02139, USA 
^ John A. Paulson Sehool of Engineering and Applied Seienees, 

Harvard University, Cambridge, Massaehusetts 02138, USA 
(Dated: March 21, 2016) 

We use scanning tunneling microscopy (STM) and quasiparticle interference (QPI) imaging to 
investigate the low-energy orbital texture of single-layer FeSe/SrTiOs. We develop a T-matrix model 
of multi-orbital QPI to disentangle scattering intensities from Fe 3dxz and 3dyz bands, enabling the 
use of STM as a nanoscale detection tool of nematicity. By sampling multiple spatial regions of a 
single-layer FeSe/SrTiOs film, we quantitatively exclude static xzfyz orbital ordering with domain 
size larger than = 20 nm x 20 nm, xzfyz Fermi wave vector difference larger than dk — 0.014 
TT, and energy splitting larger than SE = 3.5 meV. The lack of detectable ordering pinned around 
defects places qualitative constraints on models of fluctuating nematicity. 


I. INTRODUCTION 

Since the 2012 discovery of enhanced high-temperature 
superconductivity in single-layer FeSe/SrTiOs [1], the 
quest to reproduce, understand, and extend this find¬ 
ing remains urgent. Single-layer FeSe weakly coupled 
to bilayer graphene is non-superconducting down to 2.2 
K [2], but when deposited on SrTiO3(001), exhibits a su¬ 
perconducting transition temperature Tc up to 65 K [3-6] 
or 109 K [7]. Efforts to elucidate the microscopic mech¬ 
anisms behind this transformation have presently led to 
divergent viewpoints [8, 9]. At the crux of the debate 
is whether single-layer EeSe/SrTiOs exemplifies a novel 
pairing mechanism involving cross-interface phonon cou¬ 
pling, or whether it shares a common electronic mecha¬ 
nism with other iron chalcogenides already seen. 

Indications of the first viewpoint were brought forth 
by angle-resolved photoemission spectroscopy (ARPES) 
measurements, which revealed that the primary elec¬ 
tronic bands possess faint “shake-off” bands offset by 
100 meV [5, 10]. The replication of primary band fea¬ 
tures without momentum offset suggests an electron- 
boson coupling sharply peaked at qr ^ 0. The boson was 
initially hypothesized to be an O phonon mode and sub¬ 
sequently observed on bare SrTiO3(001) [11]. Model cal¬ 
culations have demonstrated that phonons can enhance 
spin-fluctuation-mediated pairing in EeSe [5, 12, 13]. 
Others have argued that phonons alone can account for 
a significant portion of the high Tc [14, 15]. 

An alternative but possibly complementary viewpoint 
is that electron doping underlies the primary enhance¬ 
ment of Tc in single-layer EeSe/SrTi03. Early exper¬ 
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iments observed that as-grown films become supercon¬ 
ducting only after a vacuum annealing process [3]. This 
procedure presumably generates interfacial O vacancies 
donating electron carriers [16]. More recent experiments 
showed that multilayer EeSe, which does not exhibit 
replica bands from coupling to SrTi03 phonons, can still 
develop superconductivity (Tc up to 48 K suggested by 
ARPES) when coated with K atoms [17-22]. Two obser¬ 
vations from these latter experiments are crucial. Eirst, 
the dome-shaped evolution of Tc with doping refocuses 
attention on electronic (spin/orbital) mechanisms of pair¬ 
ing [9]. Second, the enhanced Tc emerges from a parent, 
bulk nematic phase, characterized in multilayer EeSe by 
a small orthorhombic distortion [23] and a large splitting 
of the Ee 3dxz and 3dyz bands [24-28] . 

Nematic order, defined more generally as broken rota¬ 
tional symmetry with preserved translational symmetry, 
is a hallmark of the parent phase of iron-based super¬ 
conductors. Importantly, both spin and orbital fluctua¬ 
tions that are candidate pairing glues can condense into 
parent nematic order [29]. Eurthermore, q ^ 0 nematic 
fluctuations that extend beyond phase boundaries can 
enhance Tc [30, 31]. This mechanism operates in any 
pairing channel, with increased effectiveness in a 2D sys¬ 
tem [32] . Recent DET calculations have shown that bulk 
and single-layer EeSe exhibit a propensity towards shear¬ 
ing, but that strong binding to cubic SrTi03 suppresses 
this lattice instability [14]. It is tempting to ask whether 
in addition to suppressing nematic order, this binding 
may push the heterostructure closer to a nematic quan¬ 
tum critical point, with intensified fluctuations. 

To investigate the possible role of nematicity, we use 
scanning tunneling microscopy (STM) and quasiparti¬ 
cle interference (QPI) imaging. By generating scatter¬ 
ing through moderate disorder, quasiparticle attributes 
such as spin/orbital/valley texture, or the supercon¬ 
ducting order parameter, are manifested in selection 
rules that underlie the interference patterns. STM also 
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affords dual real- and momentum-space visualization 
of electronic states within nanoscale regions. Previ¬ 
ous STM works have uncovered C2 electronic patterns 
in parent Ca(Fei_a^Coa^) 2 As 2 [33, 34], LaOFeAs [35], 
NaFeAs [36, 37], and superconducting orthorhombic 
FeSe [38, 39]. In addition, remnant nematic signa¬ 
tures were detected in the nominally tetragonal phases 
of NaFeAs [36] and FeSeo. 4 Teo .6 [4:0]. These latter ob¬ 
servations motivate our present investigation. Can lo- 
eal disorder or anisotropie perturbations pin nanoseale 
patehes of otherwise-fluetuating nematieity in single¬ 
layer FeSe/SrTiOs, signaling proximate nematie quan¬ 
tum eritieality? Or is the heterostrueture too far from 
a nematie phase boundary for fluetuations to persist and 
boost Tc? 

This paper is organized as follows: Section II presents 
experimental details, including QPI images acquired on 
single-layer FeSe/SrTiOs. In order to extract the low- 
energy orbital texture and disentangle scattering inten¬ 
sities involving Fe Mxz and ‘4dyz bands, we develop a 
T-matrix model of multi-orbital QPI, with results shown 
in Sec. Ill and mathematical details given in Appendix A. 
In Sec. IV, we sample multiple spatial regions of our film, 
and based on our orbital-resolved QPI model, exclude 
static nematieity in the form of xzjyz orbital ordering. 
Within domains of size = 20 nm x 20 nm, we place 
quantitative bounds on xzjyz Fermi wave vector differ¬ 
ence {dk < 0.014 tt) and xzjyz pocket splitting energy 
(dE < 3.5 meV). The lack of detectable ordering pinned 
around impurities places qualitative constraints on mod¬ 
els of fluctuating nematieity. A discussion and summary 
of results is given in Sec. V. Additional details on local 
defect structure and fitting procedures are presented in 
Appendices B and C. 


II. EXPERIMENT 

Films of single-layer FeSe were grown epitaxially on 
SrTiO3(001) following procedures outlined in Ref. [41], 
then imaged in a home-built STM at 4.3 K. Typical su¬ 
perconducting gaps observed at this temperature were 
^14 meV [41]. Figs. l(a)-(c) present three atomically 
resolved topographies of the same area, acquired with 
different energy set points. Each bright spot corresponds 
to a surface Se atom; there are no in-plane defects in 
this region. Our images reveal that even pristine single¬ 
layer FeSe/SrTiOs displays appreciable electronic inho¬ 
mogeneity, in strong contrast to multilayer films grown 
on bilayer graphene [2, 38]. The electronic inhomogene¬ 
ity in FeSe/SrTiOs underscores the need for nanoscale 
measurements of electronic structure. 

To image QPI and extract local orbital information, 
we acquired conductance maps g{r^uj) = dljdV{r^eV) 
over regions of the film with in-plane defects (exempli¬ 
fied in Figs. l(d)-l(f)). A brief commentary on the de¬ 
fect structures is given in Appendix B. Figure 2(a) shows 
the Fourier transform amplitude \g{q,oo)\ of a map with 



FIG. 1. (color online) (a)-(c) Pristine region of single-layer 
FeSe/SrTiOa. Topographies of the same area are acquired 
with three set points, revealing background electronic disor¬ 
der: (a) —250 mV, 1.25 nA, (b) 50 mV, 250 pA, (c) 250 mV, 
1 nA. (d)-(f) Defect region of single-layer FeSe/SrTiOa. To¬ 
pography and differential tunneling conductance maps of the 
same area, revealing quasiparticle interference, (d) 100 mV, 
5 pA, (e) —20 mV, 200 pA, bias oscillation V-ms = 1.4 mV, 
(f) 20 mV, 200 pA, Vms = 1.4 mV. 


cj = 20 meV. Ring-shaped intensites appear around q = 
0, (±7r, Ttt), (0, ±27r), and (±27r, 0) due to scattering 
of Fermi electron pocket states. Previous works utilized 
ring size dispersion to map filled- and empty-state band 
structure [41], or energy- and magnetic field-dependent 
ring intensities to infer pairing symmetry from coher¬ 
ence factor arguments [43]. Here, we will examine ring 
anisotropy associated with the high-gr scattering channels 
(red and blue boxes in Fig. 2(a)). We will demonstrate 
that (1) although all the QPI rings are derived from scat¬ 
tering within and between the same electron pockets, the 
high-gr scattering channels have more stringent selection 
rules and hence a cleaner orbital interpretation; (2) the 



FIG. 2. (color online) (a) Fourier transform amplitude 
\g(q,(ju)\ of a 20 nm x 20 nm conductance map. Set point: 
20 mV, 200 pA; bias oscillation Vms = 1.4 mV. Drift- 
correction [42] and mirror symmetrization along the Fe-Fe 
axes are applied for increased signal. Note that fourfold ro¬ 
tational symmetrization has not been applied, (b) Grystal 
structure of single-layer FeSe. The black/gray Se atoms la¬ 
beled +/— lie above/below the plane. Dashed lines in (a) and 
(b) refer to the 1-Fe UG. 
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high-gr scattering channels can be utilized to search for 
signatures of xzjyz orbital ordering. 

Directly from the data in Fig. 2(a), we observe an un¬ 
usual relationship between the anisotropic rings around 
q = (0, 27r) and (27r, 0). In a single layer of FeSe, the 
Fe atoms are arranged in a planar square lattice, from 
which we define a 1-Fe unit cell (UC) [Fig. 2(b)] and 
crystal momentum transfer q [Fig. 2(a)]. We emphasize 
the distinction between the /c-space Brillouin zone and 
the gr-space crystal momentum transfer that is directly 
detected by STM imaging of QPI patterns. Including 
the Se atoms staggered above and below the Fe plane, 
the primitive UC becomes doubled. We might expect the 
QPI rings around (0, 27r), (27r, 0) to be identical trans¬ 
lations by the 2-Fe UC reciprocal lattice vector 2G = 
(—27r, 27r). Instead, they appear to be inequivalent and 
related by 90° rotation. The cause and implications of 
this observation will be discussed in the following section. 


III. MULTI-ORBITAL QUASIPARTICLE 
INTERFERENCE 

We develop a model to map experimental \g{q^ oo) \ pat¬ 
terns to the orbital characters of the scattered quasipar¬ 
ticles, similar in concept to Ref. [44]. In this section, we 
present an intuitive picture, followed by T-matrix simu¬ 
lations with and without xzjyz orbital ordering. Model 
details are given in Appendix A. 

Since the Fermi surface (FS) of single-layer FeSe is 
derived from Fe Sd orbitals, a natural starting point is 
to consider a low-energy model of a square lattice of Fe 
atoms. Figure 3(a) depicts a schematic FS, consisting of 
single elliptical electron pockets around k = (0, tt), (tt, 
0). The hole pockets which typically appear around (0, 
0) in other iron-based superconductors are sunken below 
the Fermi energy due to electron doping from SrTiOs [3- 
5, 45]. 

Although the Se atoms positioned between next- 
nearest neighbor Fe atoms contribute little spectral 
weight to the FS, their presence alters crystal symme¬ 
try and cannot be ignored. Their staggered arrangement 
doubles the primitive UC, folding the electron pockets 
around (0, tt), (tt, 0) on top of each other [Fig. 3(b)]. 
However, an underappreciated fact is that the pockets 
do not become identical replicas. Fe 3d orbitals that are 
even with respect to ^-reflection (x^— xy^ 3z^—r‘^) 
cannot distinguish whether Se atoms lie above/below the 
plane; only odd orbitals (xz, yz) feel an effective poten¬ 
tial of doubled periodicity [46] . In terms of tight-binding 
(TB) models, the only hopping terms that get folded in 
A;-space are those involving a product of odd and even 
orbitals [47-49]. 

To illustrate, we simulate the FS for single-layer 
FeSe/SrTiOs and show the dominant orbital contribu¬ 
tions in Figs. 3(c)-3(f). Due to incomplete folding, the 
orbital textures of the pockets around (0, tt) and (tt, 0) 
remain distinct and separately dominated by xzjxy and 




FIG. 3. (color online) (a) Schematic Fermi surface (FS) of 
single-layer FeSe/SrTiOs, derived from Fe 3d orbitals. The 
FS is composed of single elliptical electron pockets around k 
= (0, tt), (tt, 0). (b) Upon introducing the potential of stag¬ 
gered Se atoms, the electron pockets would naively fold onto 
each other (double-headed arrow); however, due to incomplete 
folding, they remain orbitally distinct, (c) Simulated FS of 
single-layer FeSe/SrTiOa and (d)-(f) its dominant orbital con¬ 
tributions {xz, yz, xy). The arrows mark the expected elastic 
scattering wave vectors which may contribute to the g = (0, 
27r) and (27r, 0) scattering channels of interest. 


yzjxy quasiparticles respectively. 

In the presence of disorder, elastic scattering channels 
should peak around wave vectors q connecting FS seg¬ 
ments with large density of states. Considering only 
the pockets shown in Figs. 3(d)-3(f), we anticipate the 
xz quasiparticles to scatter predominantly around q = 
(0, 27r), the yz quasiparticles to scatter predominantly 
around q = (27r, 0), and the xy quasiparticles to scat¬ 
ter around both wave vectors. Figures 4(a)-4(c) show 
T-matrix calculations of the orbital-resolved, density-of- 
states (DOS) modulations \pmm{Qi^ = 0)|- The index 
m denotes the Fe 3d orbitals, and we assume a local¬ 
ized, 5-wave scatterer in our simulations. Comparing 
simulation results to experimental QPI patterns at low 
energies [Fig. 2(a)], we observe that the elliptical rings 
around q = (0, 27r), (27r, 0) resemble the xz- and yz- 
projected DOS modulations respectively. Signatures of 
the x^-projected DOS modulations, which involve oppo¬ 
sitely oriented elliptical rings [Fig. 4(c)], appear to be 
suppressed in Fig. 2(a). Due to the in-plane orientation 
of xy orbitals, their wave function amplitudes at the STM 
tip height are likely smaller. 

The disentangling of derived QPI signals around 

q = (0, 27r), (27r, 0) and the suppression of xy signals yield 
a straightforward prescription to detect nanoscale xzjyz 
orbital ordering. (In contrast, the QPI signal around 
q = (0, 0) would involve both xz and yz orbital con¬ 
tributions). Such orbital ordering would lead to a pop¬ 
ulation imbalance of xzjyz carriers, implying unequal 
Fermi pocket sizes and resulting anisotropy between the 
(0, 27r), (27r, 0) scattering channels. To simulate this ef¬ 
fect, we add on-site ferro-orbital ordering to our model. 
























4 


(a)__- i-A. 

1 ^ T XZ 

(b)... 

i 

\yz 

1 xy 



4 . 

Lx 


1 

4 - 




1 

L___ 4 I| _ 





(e). 




1 


^ n • 

L -O- 


FIG. 4. (color online) Quasiparticle interference simulations 
of single-layer FeSe/SrTiOa. (a)-(c) Orbital-resolved, density- 
of-states modulations \pmm{Q^^ — 0)| in the presence of a 
localized, s-wave scatterer. The boxes mark signals around 
q = (0, 27r), (27r, 0) that are the focus of this paper, (d)- 
(f) = 0)1 with the inclusion of xzjyz orbital order. 

The arrows mark the resulting distortion of the rings. 


Figures 4(d)-4(f) show simulation results, which demon¬ 
strate squishing of the xz ring signal around (0, 27r) and 
the rounding of the yz ring signal around (27r, 0). We 
show the x^-projected DOS modulations for complete¬ 
ness, although its associated tunneling amplitude is sup¬ 
pressed. Recent studies have also proposed orbital order¬ 
ing to be bond-centered and d-wave [27, 50], but these 
complexities produce the same qualitative effect for QPI 
involving the electron pockets only. 


IV. EXPERIMENTAL BOUNDS ON 
NANOSCALE ORBITAL ORDERING 

We carried out experimental tests for xzjyz orbital 
ordering as follows. To account for local inhomogeneity, 
we sampled QPI over four distinct domains of size = 
20 nm X 20 nm (called Areas A through D in Fig. 5). 
Each domain was imaged following a separate STM tip- 
sample approach, and was likely separated from other 
domains by distances larger than our scan frame width, 
8L = 1.5 yum. To rule out tip anisotropy artifacts, the 
data from each domain were acquired with a different 
microscopic tip termination, modified by field emission 
on polycrystalline Au. Over every domain, conductance 
maps were acquired at low energies (±10 meV) in order 
to compare scattering from xz and yz Fermi pockets. 

Figures 5(a) and 5(b) show QPI images acquired over 
area A, cropped around q = (0, 27r) and (27r, 0). We ap¬ 
plied Gaussian smoothing with width a = where ^q = 
0.028 TT is the inherent resolution for momentum defined 
within a finite 20 nm x 20 nm area. To compare the xz- 
and ^ 2 ;-derived QPI rings, we took line cuts along their 
minor axes (arrows in Figs. 5(a) and 5(b)), where the sig¬ 
nal intensity is the strongest. Each line cut was averaged 
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FIG. 5. (color online) Nanoscale wave vector analysis of or¬ 
bital ordering, (a), (b) Gropped quasiparticle interference 
images \g{q,uj)\ around (0, 27r) and (27r, 0), with arrows in¬ 
dicating line cuts used to compare xzjyz Fermi pocket sizes 
(perpendicular bars represent averaging width). (c), (d) Line 
cuts from conductance maps acquired over four distinct do¬ 
mains (8r^ — 20 nm x 20 nm), labeled A through D, and 
two energies, ±10 meV. The horizontal bars mark the peak 
locations determined from Gaussian fits (solid lines), with in¬ 
herent resolution 8q — 0.028 tt. For visualization, the line 
cuts are normalized by the Bragg peak amplitude and verti¬ 
cally offset. 


over a width of 10 pixels. The results for the four do¬ 
mains and two energies are shown in Figs. 5(c) and 5(d). 
We determined peak locations from Gaussian fits (solid 
lines). The horizontal bars denote extracted peak loca¬ 
tions with error ±(5g. The addition of a linear background 
is found to shift the fitted Gaussian peak locations, by 
an amount smaller than 8q (Appendix C). In all four ar¬ 
eas, we observe no significant deviations between the xz- 
and ^z-derived QPI wave vectors. We therefore exclude 
orbital ordering with domain size larger than = 20 
nm X 20 nm and xzjyz FS wave vector difference larger 
than = (5g/2 = 0.014 tt (the factor of two arises when 
changing between q space and k space). 


We also determine an energy bound on xzjyz orbital 
ordering. Figures 6(a) and 6(b) show a simulated split¬ 
ting of the xzjyz bands for reference, and Figs. 6(c) and 
6(d) show the corresponding QPI dispersions measured 
over Area B. The dispersing peaks locations are extracted 
from Gaussian fits, shown in Fig. 6(e), and are identical 
within ±(5g over the given energy range [—30 meV, 20 
meV]. Due to the overlap with a sunken, zone center 
hole pocket [41, 45], the lower edges of the Fermi pock¬ 
ets are difficult to detect. Instead, we fit the dispers¬ 
ing peak locations to parabolas, and find their respective 
band edges to be —51.5±3.5 meV and —49.6±3.0 meV. 
We again bound orbital ordering with lyz < 5E = 
3.5 meV. 
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FIG. 6. (color online) Nanoscale dispersion analysis of orbital 
ordering, (a), (b) Simulated pockets with xzjyz orbital order¬ 
ing, revealing a split in the band edges of magnitude 2/\xziyz 
— 60 meV. (c), (d) Experimental intensity plots of quasipar¬ 
ticle interference images \g{q^uj)\ acquired over Area B, cut 
along the arrows in Figs. 5(a) and 5(b). (e) Plot of dispersing 
peaks positions extracted from Gaussian fits. The width of 
the horizontal bars is = ±0.028 tt. The superconducting 
gap magnitude is marked by 2A. 


V. DISCUSSION AND SUMMARY 

We return to the central debate, whether single-layer 
FeSe/SrTiOs exemplifies a novel interface-phonon pair¬ 
ing mechanism, or whether it can be explained by 
an electronic mechanism common to other iron chalco- 
genides. Recent experiments demonstrating sizeable Tc 
in electron-doped multilayer FeSe [17-22] suggest an elec¬ 
tronic pairing mechanism in other iron chalcogenides 
lacking hole Fermi surfaces, distinct from a SrTiOs 
phonon mode. One possibility involves qr ^ 0 nematic 
fluctuations extending from the parent ordered phase. 
Theories have shown that such fluctuations can boost Tc 
effectively in any pairing channel, on both the ordered 
and disordered sides of the phase transition [30-32]. 

Moving from multilayer to single-layer FeSe/SrTiOs, 
we face two scenarios. One scenario is that single-layer 
FeSe/SrTiOs remains close to a nematic phase boundary. 
Here, ordering is absent, but intense nematic fluctuations 
may be pinned by impurities. Another scenario is that 
single-layer FeSe/SrTiOs lies sufficiently far away from 
a nematic phase boundary, such that C2 electronic sig¬ 
natures are not produced even upon local perturbation. 
The quantitative bounds on static xzjyz orbital ordering 
derived from our QPI measurements favor the latter sce¬ 
nario. In turn, this statement would suggest that single¬ 
layer FeSe/SrTiOs is not the same as electron-doped mul¬ 
tilayer FeSe, in which nematic fluctuations may be oper¬ 
ative [19]. The addition of the SrTiOs substrate intro¬ 
duces novel effects beyond electron doping, such as pos¬ 
sible interface phonons, that push the two systems apart 
in phase space. 

To summarize, we have utilized STM and QPI imag¬ 
ing to demonstrate that the pronounced nematic order 
present in multilayer FeSe is suppressed in single-layer 


FeSe/SrTiOs. More importantly, nanoscale nematic or¬ 
dering is not recovered upon perturbation by anistropic 
defects. We arrived at our conclusions by comparing 
high-qr scattering channels around (0, 27r) and (27r, 0), 
which we showed by T-matrix simulations to be sepa¬ 
rately dominated by xz and yz quasiparticles. Our work 
places quantitative bounds on static xzjyz orbital or¬ 
dering in single-layer FeSe/SrTiOs, and qualitative con¬ 
straints on models of Tc enhancement by nematic fluctu¬ 
ations. 
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Appendix A: Model of multi-orbital quasiparticle 
interference 

Model Hamiltonian: We begin with a TB model for 
single-layer FeSe where the low-energy bands are pro¬ 
jected onto the five 3d orbitals of an Fe atom: 

^ imn{\'lx ~ jx\^ \'ly ~ I )^m (j ) • (-^1) 

ij mn 

Here, i, j index the Fe lattice sites and m, n index the five 
orbitals. The tilde symbol indicates that a momentum 
shift Q = (tt, tt) has been applied to the even orbitals in 
order to downfold the UC from two Fe atoms to one [48, 
49, 51, 52]. The corresponding bare Green’s function is 
given by 

&{k,i0)= [{io + i8)l5y,5 - (A2) 

where the bolded capital symbols are matrices and 6 is a 
broadening (= 5 meV for all simulations). 

We adapt hopping parameters tmn computed in 
Ref. [53] , then apply rescaling to qualitatively capture the 
low-energy spectrum of single-layer FeSe/SrTiOs [41]. 
The hopping terms are given in Table I, and the resulting 
band structure is shown in Fig. 7. The electron pocket 
around k = (0, tt) remains attached to a hole pocket. 
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mn = 55 

-0.1504 
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-0.0058 


0.0135 

mn = 13 
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0.0381 


0.0081 



mn = 14 


0.1304 


0.0054 


0.0108 



mn = 15 


-0.0762 


-0.0327 




-0.0054 

mn = 34 






-0.0038 



mn = 35 


-0.1154 




-0.0077 



mn = 45 




-0.0577 




0.0038 


TABLE L Rescaled hopping parameters for tight-binding model adapted from Ref. [53]. Here, m—1 is xz^ m—2 is yz^ m=3 is 

22 A * r'*o22 

X -y , m—A IS xy^ m—o is 6z -r . 


but this does not affect our simulation results closer to 
the Fermi level. The positions of the T pockets above and 
below the Fermi energy also do not affect our simulations. 

Fermi surface simulations: Computing the FS of 
single-layer FeSe/SrTiOs [Figs. 3(c)-3(f)] from our TB 
model requires that we restore the original crystal sym¬ 
metry (with a 2-Fe UC) induced by staggered Se atoms. 
To do so, we transform the lattice operators with a site- 
dependent sign factor [48]: 

cl{i) = (A3) 

Here, Pm = ±1 for orbitals that are even/odd with re¬ 
spect to 2 :-reflection. This transformation is equivalent to 
undoing the downfolding operation applied in Eq. A1 and 
shifting the even orbitals by —Q in momentum space: 


4n{^) = 


^m(^) 5 


Pm — I5 
Pm ~ 


(A4) 


where kq = k — Q. The orbital components of the FS 



FIG. 7. Band structure of single-layer FeSe/SrTiOa in the 1- 
Fe Brillouin zone. The tilde symbol indicates that a momen¬ 
tum shift Q = (tt, tt) has been applied to the even orbitals 
to downfold the unit cell from two Fe atoms to one. Adapted 
from Ref. [53], with rescaled hopping parameters to match 
the pocket edges measured in Ref. [41]. 


are given by 


= 0) 


Pm — 

Pm = +!> 


(A5) 


where A^^(k^uj) = —Im G^^{k^uj)/7r. Further insights 
on (1) the connection between Eq. A5 and ARPES- 
measured band structures, (2) common misconceptions 
of whether spectroscopic probes measure quasiparticles 
closer to the 1-Fe or 2-Fe Brillouin zone description, and 
(3) proper folding of the superconducting pairing struc¬ 
ture, are given in Refs. [47-49]. 

Quasiparticle intereference simulations: To generate 
QPI, we introduce a localized, s-wave scatterer at i = (0, 
0) of uniform strength V = 1 eV in all orbital channels. 
The resulting impurity Green’s function is given by 

G{k,k',uj) = G\k,uj)T{uj)G\k',uj), (A6) 


for k ^ k'^ and the T-matrix is momentum-independent: 


T{w) = 


-^5x5 



cPk 


G°{k,Lu) 


-1 


V. 


(A7) 


Since STM measures local density of states in real 
space, we additionally transform lattice operators cj„(i) 
into continuum operators 


V’m(^) = “ *)cL(*)- (AS) 

i 


The first factor on the right recovers the proper crystal 
symmetry (2-Fe UC) due to staggered Se atoms [Eq. A3]. 
The second factor on the right, 0^, is the Wannier func¬ 
tion associated with orbital m at site i. This factor 
captures nonlocal tunneling contributions [54]. For sim¬ 
plicity, we approximate the Wannier functions at the 
STM tip height with a square cutoff in momentum space: 
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Uk)AO^Jk, Oj \p^Jq, 0;| 



FIG. 8. (color online) (a) Lattice model, in which the Green’s 
function has nonzero weight restricted to discrete lattice 
points; i.e., the Wannier functions are given by 0(r — i) — 
S{r — i). (b) In momentum space, (l)m{k) = 1, such that 

there is no cutoff for states involved in scattering. Shown 
here is (l)rn{k)A^^{k,uj = 0) for the xz orbital, (c) Gonse- 
quently, additional ring intensities arise in \pmm{q, 0)| around 
q = (27r, 0), (27r, 27r) (arrows in (b), boxes in (c)) that are not 
observed experimentally, (d) Gontinuum model, which in¬ 
corporates nonlocal effects due to a finite Wannier function 
width (white square), (e) We model the experimental data 
with Wannier functions of the form 0m (fc) = 1 for kx, ky G 
[-I.Stt, I.Stt], and <j)m(k) — 0 otherwise (white square), (f) As 
a result, there are fewer scattering channels. 


0m(^) = 1 foi* ky G [-I.Stt, I.Stt], and 0m(^) = 0 oth¬ 
erwise. In real space, this corresponds to a characteristic 
tunneling width of 0.67 (ape-Fe), which is neeeded to re¬ 
produce experimental QPI patterns. Figure 8 illustrates 
qualitative differences between simulations with (contin¬ 
uum model) and without (lattice model) non-local tun¬ 
neling. 

From Eq. A8, we obtain the continuum impurity 
Green’s function: 

Qmm (^7 k ^ 

_ i k , a;)0^(/c)0m(^ )? Pm — (^Q) 

^ Gmm ), Pm — + 1 . 

Only diagonal elements are shown for brevity. Finally, 
the orbital projections of the DOS modulations are given 
by 


Pmm{qi — r, 


/ 


dPk 


2tt J (27r)^ 


0mm (^7 k q^Cj) 

- Qmmik,k - q,u) 


(AlO) 


Figs. 4(a)-(c) show plots of |pmm (^7 0)| for the xz^ yz, 
and xy orbitals. 

Orbital ordering: To simulate on-site, ferro-orbital or¬ 


dering, we include the following term [55] in our TB 
Hamiltonian [Eq. Al]: 

^xzjyz ~ ^xzjyz (^) (^)Cy2; (^)] • (^fl) 

i 

A value of A^^jyz = 30 meV was used for Eigs. 4(d)-(f) 
and 6(a)-(b). 

Superconductivity: The inclusion of superconductiv¬ 
ity does not change the QPI orbital texture. Eollow- 
ing Ref. [56], we introduce superconductivity in band 
space, but compute scattering in orbital space. Erom the 
normal-state TB Hamiltonian [Eq. Al], we define bands 
e{k) = U{k)H^{k)U^(k), where U{k) represents a uni¬ 
tary transformation. The Green’s function in the super¬ 
conducting state is then given by 

G°cs(fe,^)= [(^ + *<5)Foxio-^^cs(fe)]“\ (A12) 

where H^Q^{k) has the following form in Nambu repre¬ 
sentation: 

^BCs(^) = 

f U^k)e{k)U{k) U^k)A{k)U^{-k) \ 

\U^{-k)A%k)U{k) -U^{-k)i{-k)U%-k)) ‘ 

(A13) 

We model isotropic gaps in band space: A{k) = A/sxs, 
with A = 14 meV based on our dl/dV measurements. 
(Recent ARPES measurements have detected small gap 
anisotropy [57]). We also take the impurity potential of 
a localized, non-magnetic, s-wave scatterer: 

V = 0 ) . (A14) 

f 0 -y/sxs) 

Eigure 9 shows QPI simulations with the inclusion of 
superconductivity, at energy cj = 10 meV. There is little 
difference compared with the normal-state calculations, 
without or with xzjyz orbital ordering. 

Anisotropic scatterer: In an angular momentum ex¬ 
pansion of the T-matrix, the leading component should 
be s-wave; i.e., intraorbital scattering, with Vmn = 
Vmm^mn^ should dominate. Vmm can in general vary 
with orbital, but this simply modifies the relative weights 
of the orbital-resolved DOS modulations. In Eig. 10, 
we illustrate this effect in the case of a C2 scatterer 
+ yyz,yz)> Tuning the strengths of V^^z^xz and 
yyz,yz tunes the intensity of the q = (0, 27r) and (27r, 0) 
scattering channels respectively; however, the scattering 
wave vectors remain unchanged and are a more robust 
measure of orbital ordering. On the other hand, the q 
^ (0, 0) channel will display anisotropies related to the 
scattering potential, so we do not analyze it. 

Equivalent ten-orbital formulation: We derive an 
equivalent formulation of the continuum impurity 
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FIG. 9. (color online) Quasiparticle interference simula¬ 
tions in the superconducting state, with isotropic gaps of 
14 meV on all bands, (a)-(c) Orbital-resolved spectral func¬ 
tion ^Bcs,mm(<?? 10 mcV). (d)-(f) Orbital-resolved density- 
of-states modulations |pBCS,mm(^, 10 meV)| in the presence 
of a localized, non-magnetic, s-wave scatterer. (g)-(i) Same 
as (d)-(f), but including xzjyz orbital ordering {A^ziyz — 30 
meV). 


Green’s function for ten-orbital TB models, such as those 
in Refs. [55, 58]. The ten 3d orbitals come from the 
two inequivalent Fe atoms (^4 and B) of the primitive 
UC; (xy)^, {x^ - (ixz)^, (iyz)^, (z^)^, (xy)^, 

(a;2 {-ixz)^, {-iyz)^, {z'^)^. 



FIG. 11. Quasiparticle interference simulations derived from 
calculating the joint density of states separately for each or¬ 
bital. 


orbital TB model: 

Qmnir,r',uj) = 

4>*m{r - i)4>n{r' - j)Gmn{i,3,^)- (A15) 

Here, is the lattice impurity Green’s func¬ 

tion, given in momentum space by Eq. A6. Next, we 
decompose the Fe lattice into two sublattices: 


A = {i : ix + iy = odd}, 
B = {i : ix + iy = even}. 


(A16) 


For diagonal terms {m = n), the sum in Eq. A15 can be 
split into four contributions: 


i,j ieA,jeA 

ieAjeB ieBjeA 

+ ^ (...)(+!). (A17) 

ieBjeB 

Taking the Eourier transform of Eq. A17 yields 


We begin with the real-space representation of the con¬ 
tinuum impurity Green’s function [Eq. A9] for a five- 



'4zxz='l0eV '/xzxz = 0-4eV 

Vy,y,= 0Aey Vy,y,= ^.0ey 


FIG. 10. QPI simulations in the presence of a localized, 
anisotropic scatterer. While the relative intensities of the 
xzjyz scattering channels around q — (0, 2ti) and (27r, 0) 
are modified, their wave vectors remain unchanged. 


T Prn^rn-\-b,rni}^'i ^ 7 T (^7 ^ 7 

T ^m-|-5,m-|-5 (fe,fe',w)]C(fe)0„(feO. (A18) 

Here, m runs from 1 through 5, Qmm{k^k\uj) is de¬ 
rived from a five-orbital TB model, and Gmn(^ 7 ^^^) is 
the lattice Green’s function for a ten-orbital TB model. 
The Pm factors appear because of minus signs present in 
the orbital definitions of {—ixz)^^ {—iyz)^. The middle 
terms in Eq. A18, which mix orbitals m and m + 5, rep¬ 
resent intraorbital basis site interference. Importantly, 
the sum of these terms are non-zero for a finite Wan- 
nier function width. These crucial terms, which have not 
been considered in previous ten-orbital QPI calculations 
of Ee-SCs [56], are required in order to reconcile five- 
orbital and ten-orbital QPI calculations in the presence 
of non-local tunneling. 

Joint density of states: Eigure 11 demonstrates that 
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FIG. 12. (color online) Differential tunneling conductance 
map revealing dominant type of anisotropic impurities (re¬ 
produced from Fig. 1(e)). The defects are directed along the 
crystalline Fe-Se axes and appear in four possible orientations 
(yellow arrows). Set point: —20 mV, 200 pA, bias oscillation 
Vrms = 1.4 mV. 

our T-matrix simulations can be qualitatively approxi¬ 
mated by calculating the joint DOS separately for each 
orbital: 

j w)A%,^{k + q, uj)4>*^{k)4>m{k + q). 

(A19) 

Appendix B: Symmetry breaking in local defect 
structures 

In NaFeAs [36, 37], multilayer FeSe [38], and 
Ca(Feia,Coa,) 2 As 2 [34], the dominant atomic-scale defects 
have been shown to pin larger electronic dimers that 
are unidirectional within nanoscale domains and aligned 
along one Fe-Fe axis, providing evidence of local nematic 
ordering. Here we search for a similar effect in single¬ 
layer FeSe/SrTiOs. 

A closer inspection of Fig. 1(e) reveals a dominant in¬ 
plane defect that appears as adjacent bright and dark 
atoms along the crystalline Fe-Se axes [Fig. 12], and ex¬ 
ists along all four orientations, like impurities observed in 
LiFeAs [59, 60]. These defects are similar in their atomic- 
scale structure and Fe-Se orientation to the dominant de¬ 
fects in NaFeAs [36, 37] and multilayer FeSe [38]. How¬ 
ever, the defects observed in single-layer FeSe/SrTiOs 
do not show the larger Fe-Fe electronic dimers. Further¬ 


more, in Sec. IV we considered four different nanoscale 
domains (20 nm x 20 nm), each containing several such 
randomly-oriented defects, but our nanoscale wave vector 
and dispersion analyses found no significant difference in 
electronic response between the Fe-Fe axes. The chance 
that the impurity orientations would exactly balance in 
all four sampled regions is small. 

Appendix C: Fitting details 

Here we detail the fitting procedures used to derive 
wave vector bounds on nanoscale orbital ordering [Figs. 5 
and 6]. We performed Gaussian fits using the iterative 
Levenberg-Marquardt algorithm implemented by Igor 
Pro. Fit errors were estimated from residuals and rep¬ 
resent one standard deviation of the fit coefficient from 
its true value, assuming independent and identically- 
distributed Gaussian noise. As seen in Table H, the fit 
errors {Sqo) for go are insignificant compared to the data 
resolution Sq. Note that gives the correlation length 
of the QPI patterns, rather than the uncertainty in its 
wave vector. 

Since there is a sizeable slope in the line cuts, we also 
performed Gaussian fitting with linear background for 
comparison [Fig. 13 and Table H]. We find that the go 
values are shifted by amounts smaller than the data res¬ 
olution Sq. The leading error is therefore the data reso¬ 
lution Sq, which we report as short horizontal blue and 
red bars in Figs. 5(c) and (d). 



FIG. 13. (color online) Line cuts from conductance maps, re¬ 
produced from Figs. 5(c) -(d). The solid lines denote Gaussian 
fits with linear background. Fit parameters are recorded in 
Table II. 
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Line cut label 

Gaussian plus constant background 

Gaussian plus linear background 


Qo ± 5qo [2 tt] 

aq± Sag [27r] 

qo±SqQ [2 tt] 

aq± Sag [2 tt] 

Area A, 10 meV, xz 

0.145T0.002 

0.08±0.70 

0.1506T0.0007 

0.021T0.01 

Area A, 10 meV, yz 

0.152T0.001 

0.017±0.004 

0.1584T0.0004 

0.0238T0.0006 

Area B, 10 meV, xz 

0.146T0.001 

0.08±0.60 

0.1518T0.0006 

0.0233T0.0009 

Area B, 10 meV, yz 

0.133T0.002 

0.05±0.04 

0.139T0.002 

0.027T0.003 

Area C, 10 meV, xz 

0.1456T0.0005 

0.021±0.004 

0.148T0.003 

0.029T0.004 

Area C, 10 meV, yz 

0.1429T0.0009 

0.1±0.5 

0.152T0.001 

0.017T0.001 

Area A, —10 meV, xz 

0.125T0.001 

0.026T0.009 

0.1307T0.0008 

0.021T0.001 

Area A, —10 meV, yz 

0.130T0.002 

0.02T0.01 

0.1363T0.0004 

0.0257T0.006 

Area B, —10 meV, xz 

0.1300T0.0003 

0.023T0.002 

0.1324T0.0006 

0.0228T0.0008 

Area B, —10 meV, yz 

0.1217T0.0007 

0.033T0.007 

0.1281T0.0009 

0.022T0.001 

Area D, —10 meV, xz 

0.130T0.001 

0.018T0.003 

0.1361T0.0009 

0.023T0.001 

Area D, —10 meV, yz 

0.1275T0.0006 

0.022T0.004 

0.129T0.001 

0.025T0.002 


TABLE IL Comparison of fit parameters between (1) Gaussians with constant background, g — Qq + Aexp(—(g — /{2a‘l)) 

[Figs. 5(c)-(d)], and (2) Gaussians with linear background, g — Qq -\- mx + Aexp(—— ^o)^/(2crg)) [Fig. 13]. As reference, the 
data resolution 5q — 0.014 27r. 
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